Critical Behavior in Trapped Strongly Interacting Fermi Gases 
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We investigate the width of the Ginzburg critical region and experimental signatures of critical 
behavior in strongly interacting trapped Fermi gases close to unitarity, where the s-wave scattering 
length diverges. Despite the fact that the width of the critical region is of the order unity, evidence 
of critical behavior in the bulk thermodynamics of trapped gases is strongly suppressed by their 
inhomogeneity. The specific heat of a harmonically confined gas, for instance, is linear in the 
reduced temperature t — {T — Tc)/Tc above Tc. We also discuss the prospects of observing critical 
behavior in the local compressibility from measurements of the density profile. 
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I. INTRODUCTION 

Atomic Fermi gases close to unitarity fl*] are among the 
most strongly interacting systems known. Understanding 
their thermodynamic properties is an iniportant theoret- 
ical [l,[l,[l,[EPl and experimental 0,11, HI challenge that 
is only beginning to be addressed. In typical strongly in- 
teracting superfluids such as ^He, thermodynamic mea- 
surements at low temperatures outside the critical re- 
gion reveal the nature of the elementary quasiparticles 
(e.g., phonons and rotons) fui\. Conversely, in the criti- 
cal region close to the superfluid transition temperature 
Tc, phase fluctuations of the order parameter play a cru- 
cial role and simple theories based on this quasiparticle 
picture breakdown. In this region, the temperature de- 
pendencies of thermodynamic quantities are instead gov- 
erned by universal scaling laws The width of the 
critical region reflects the strength of interactions in the 
system: for weak-coupling BCS superconductors, it is 
very small, with magnitude {Tc/epY 10~^^ [l2]. In 
contrast, for strongly interacting clean high- Tc cuprate 
superconductors, it can reach 10~^ [H, lU- In trapped, 
weakly interacting atomic Bose gases, it is typically be- 
tween 10""* and 10-2 One of the widest crit- 
ical regions belongs to superfluid ^He, where it is of the 
order unity. The existence of such a large critical region 
gives rise to a clear signature of critical behavior in ther- 
modynamic quantities, most dramatically in the specific 
heat which exhibits the famous lambda curve p^ . 

An important question in the study of trapped strongly 
interacting Fermi gases is to what extent the measured 
temperature dependencies of thermodynamic quanti- 
ties 0, S [§] reflect a gas of weakly interacting quasi- 
particles as opposed to universal scaling behavior aris- 
ing from phase fluctuations in the critical region. In 
contrast to the situation in trapped weakly interacting 
Bose gases, the strong interactions in Fermi gases close 



to unitarity can give rise to a large critical region, sug- 
gesting that scaling behavior (of the three dimensional 
XY universality class) might be evident in bulk thermo- 
dynamic quantities such as the specific heat. Interpreting 
thermodynamic measurements in trapped gases is com- 
plicated by the inhomogeneity of the gas, however, since 
the superfluid transition temperature Tc{r) depends on 
the distance r from the center of the trapping potential 
(assumed to be isotropic). This means that sharp signa- 
tures of critical behavior in bulk thermodynamic quanti- 
ties that exist in homogeneous samples become smoothed 
out by the averaging effects of the nonuniform density 
profile in the trap. 

Our aim in this paper is to understand whether or not 
current experiments with trapped Fermi gases are ob- 
serving signatures of critical behavior and also to act as 
a guide for future experiments looking for such signa- 
tures. Building on the work of Sa de Melo, Randeria, 
and Engelbrecht ^19'|, in Sec. |TT1 we study the width of 
the critical region of a two-component Fermi gas through 
the BCS-BEC crossover [2^ using the microscopic theory 
of Nozieres and Schmitt-Rink (NSR) [2l|. That the crit- 
ical region in a Fermi gas close to unitarity should have 
a width of order unity has been pointed out by Rande- 
ria [13], based on the fact that the T = Ginzburg- 
Landau coherence length is very small there, comparable 
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to the mean inter-particle spacing fcp [23|, |2j| . We con- 
firm that the critical region has a width of order unity in 
the strongly interacting region close to unitarity, reaching 
a maximum just on the BEC side of resonance. 

In Sec. IIII[ we derive the universal scaling relation 
for the trap averaged bulk specific heat C = J drc{r) 
in a Fermi gas confined by a power-law trapping poten- 
tial. We also discuss the scaling properties of the local 
isothermal compressibility Kxir), since this quantity can 
be measured by looking at the density profile n{r) of the 
trapped gas. 

The effect of a harmonic trapping potential on bulk 
thermodynamic quantities in a Bose gas has been inves- 
tigated by Damle et al. in Ref. Using a local den- 
sity approximation, they calculate the scaling behavior 
of various quantities by substituting the local transition 



2 



temperature Tc(r) into well-known scaling relations de- 
veloped for homogeneous three dimensional XY systems. 
This approach has also been used to determine the spe- 
cific heat critical exponent a in high-precision measure- 
ments d^ll^H^I performed on ^He in large sample sizes, 
where gravity effects produce a position-dependent tran- 
sition temperature as well. In this paper, we carry out 
a detailed analysis of the bulk specific heat of a strongly 
interacting trapped Fermi gas. Despite the large criti- 
cal region, we find that the critical scaling behavior of 
the bulk specific heat in a trapped Fermi gas is strongly 
suppressed relative to that in a uniform superfluid. 

Finally, in Sec. IIVI we discuss the relevance of our 
results to current experiments on strongly-interacting 
Fermi gases. 



II. GINZBURG CRITICAL REGION 

The Ginzburg critical region defines the range of tem- 
peratures about Tc where the root mean square thermal 
fluctuations of the order parameter exceed some char- 
acteristic mean- field value [28j . In the superfluid state, 
this characteristic value is the mean-field order parame- 
ter while an analogous value can be defined in the normal 
phase. The width of the critical region is estimated from 
the Ginzburg-Landau functional ^28l] 



An= dr 



A\SA\^ + ^\SA\^ + D\VSAf 



(1) 



for fluctuations SA of the order parameter Aq. (This ex- 
pression is valid for the normal phase, where Aq = 0. A 
similar expression describes the fluctuations in the super- 
fluid phase Here, A{T) = at, where t = {T-Tc)/Tc. 
In the vicinity of Tc, the root mean square thermal fluc- 
tuations of the order parameter are much smaller than 
the characteristic value y^2at/B as long as t ^ Gi, 
where 1291 
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(2) 



is the Ginzburg-Levanyuk number. For values of T in- 
side the critical region {\t\ < Gi), fluctuations are large, 
meaning that the mean-field Landau theory of phase 
transitions breaks down and critical behavior sets in. 

In the remainder of this Section, we calculate the 
width (Gi) of the critical region through the BCS-BEC 
crossover |2l| for Fermi superfluids using NSR theory. 
For a two-component gas of fermions, tuning the s-wave 
scattering length Og, one can access the entire crossover, 
from a weak-coupling BCS superfluid of Cooper pairs 
(tts small and negative) to a Bosc-condcnsate of tightly- 
bound dimer molecules in the BEG limit {a^ small and 
positive). In between, the scattering length diverges 
to inflnity at unitarity. The region of strong interac- 
tions close to unitarity has been realized in experiments 
on trapped Fermi gases with Feshbach scattering reso- 
nances |l!|. 
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FIG. 1: (Color online) Ginzburg critical region tlirougli the 
BCS-BEC crossover. The dashed line denotes the critical 
temperature Tc calculated within the NSR formalism as a 
function of the dimensionless interaction parameter (fcpas)"^ 
from the BCS region ((fcpas)"^ < 0) into the BEC region 
((ftpfls)"^ > 0) (from the calculation in Ref. H^). The shaded 
region is the critical region defined by |T — Tc|/Tc < Gi. Inset: 
the Ginzburg-Levanyuk number Gi. 



NSR theory is the simplest theory of the BGS-BEG 
crossover that gets the correct superfluid transition tem- 
perature in the BEC limit of strongly-coupled dimer 
molecules. In contrast to mean-field BCS theory, which 
predicts a transition temperature equal to the disassocia- 
tion temperature of dimer molecules, NSR gives the tem- 
perature at which a gas of dimers Bose-condenses [l9l.[2l|. 

0.2ITf (see Fig. [T|), where Tp = er/feB is the Fermi 
temperature and ep is the Fermi energy). We emphasize, 
however, that NSR is still "mean-field" insofar as it is a 
Ginzburg-Landau expansion of the thermodynamic po- 
tential (action) in powers of order parameter fluctuations 
close to Tc: ^ 



qpp' 



/3 

(3) 

Here, q = {q^,iqm), p, and p' are momentum/Matsubara 
frequency four- vectors. 
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is the inverse two-particle vertex function for fermions 
with s-wave scattering length a^, /k = [exp(/3^k) + 1]^^ 
is the Fermi thermal distribution, ^k = Ck — ^^id ek — 
k2/2m. We set = 1 throughout this paper, unless 
stated otherwise. Uq is the long- wavelength, static limit 
of the effective interaction between Cooper pairs: 
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Here, X = tanh(/?^k/2) and Y = /?sech2(/3ek/2). Tc 
is determined from Fp^'^ly^j'^ = 0, where the chemical 
potential fx is calculated from the number equation [l^ 
n = fc|/37r2 = ^^^(1 -X)+ 13-^ T^'dT./d^i, where 
fcp is the Fermi wavevector. 

The coefficients in the classical Ginzburg-Landau 
bmctional, ([1]), can be obtained from the static limit 
{qm — 0) of the quantum action, ([3]), using AQ = S/ (3. 
They are given hy D ^ 5(9^r^V5q^)l9=o,T=T<,, a = 
Tc(9ro V^T)|t=t,, and B = Uo\t=t,. Explicitly, 
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(7) 



where /3c = 1/ fceTc and all terms are evaluated aiT = T^. 

There is no analog of the second term in ([7]) (involv- 
ing the temperature derivative of the chemical potential) 
in the classic analysis of Gork'ov [33], which considered 
the BCS limit of the BCS-BEC crossover. In the ECS 
limit, this term is negligible. However, it becomes in- 
creasingly important through the crossover, finally be- 
coming the dominant term in the BEG limit. In this 
limit, where Og > is small and positive, the NSR theory 
of a two-component Fermi gas reduces to the Bogoliubov- 
Popov theory [3l[ of a dilute Bose gas of dimer molecules 
with mass M = 2m [s^ . [ssj . In this theory, the Bose 
chemical potential vanishes at and, constructing 
a Ginzburg-Landau theory based on it, we would have 
A = -fiB and a = -TcidfiB/dT)T=T,. In the BEG limit 
of NSR theory, F^^ -> Z(-ig,„ + q^/4m - fiB ), where 
Z = m^Os/STT and ^b = 2fi + \Ei,\ is the chemical po- 
tential for dimer molecules with binding energy = 



-1/maf [311331. The Thouless criterion F^"|t=t, = 
for Fermi superfluids thus also requires fiB to vanish at Tc 
and consequently, a — —ZTc{dpLB/dT)T=T^^ consistent 
with the Bogoliubov-Popov theory of dimer molecules. 
Thus, we see how the second term in ([7]) is crucial for 
obtaining the correct behavior of a in the BEG limit. 

In the BEG limit, the chemical potential ~ i?b/2 ^ 
is approximately one-half the binding energy of a 
dimer molecule and one finds D = mas/'il-K, B = 
m^al/lQi:, and a = -[Tc{^^i/^T)\T=T^rn^as/'^^T. Ex- 
panding the BGS gap equation in powers of Ao/|/i|, one 
obtains (see also Eq. (38) in Ref. [l^]) A's — UdcindT) 
in the superfluid phase, where Udd = 'inas/m is the 
mean- field dimer-dimer interaction and [35| nc{T) = 
AQ(T)m^as/87r is the condensate density of dimcrs. Us- 
ing the approximation Ao(T) « Ao(0)-\/^ with Ao(0) = 
Aep/^/Sirkpas, one finds a = mkpa'^ /12Tr'^ and Gi ~ 
127r(fcFas)(fcBlc/eF)^ ^ {kpcis)- Note that this coincides 

1 /3 

with a known condition \t\ nj add for the Bogoliubov- 
Popov theory of a dilute Bose gas (here, a gas of dimers 
with density = n/2 and scattering length a^d — 2as) 
to break down |31i| . 



In the BGS limit where < is small 

and /i ~ ep, it is straightforward to show (see 
the related calculations in Ref. [ll|) that a = 
iV(eF), D = 7C(3)iV(eF)eF/247r2m(fcBTc)2, and B = 
7C(3)iV(eF)/87r2(fcBTc)2, where 7V(eF) is the density of 
states at the Fermi surface. Thus, the width of the criti- 
cal region is proportional to (Tc/eF)"*, as expected. 

In the region close to unitarity, the GL coefficients 
are determined numerically using the NSR values for the 
chemical potential at Tc [s^. To avoid errors associ- 
ated with taking a numerical derivative of the chemi- 
cal potential with respect to temperature (in fact, the 
temperature dependence of the NSR chemical potential 
exhibits unphysical features in the critical region [ssj). 
we take the temperature derivative of the gap equa- 
tion. This gives us an equation for {dfj./dT)\T=T^ in 
terms of ^c, Tc, and (9Ao/(3T)t=Tc in the superfluid 
phase. (Strictly speaking, the Ginzburg-Levanyuk num- 
ber will assume different values in the normal and su- 
perfluid phase. This difference is not crucial, however, 
and we take Gi to be symmetric.) For this last quan- 
tity, we use Ao(T) « Ao(0)-\/^, where Ao(0) is taken 
from Ref. [35|. We plot the critical region through the 
BGS-BEG crossover in Fig. [1] The Ginzburg-Levanyuk 
number is very small in the BGS region {{kpas)^^ < 0), 
approaches unity close to unitarity, and slowly shrinks 
again in the BEG region ((fepOs)^^ > 0). It is interest- 
ing to note that the critical region is actually widest on 
the BEG side of resonance, close to the point where the 
chemical potential vanishes. Related nonmonotonic de- 
pendence of the Ginzburg-Landau coherence length was 
found in Ref. [H. 

For a trapped gas, the fact that the Ginzburg criti- 
cal region close to unitarity is of order unity means that 
when the temperature at the center of the trap reaches 
the local superfluid transition Tc{r — 0), a significant 
fraction of the gas will be in the critical region. In the 
critical region, thermodynamic quantities should exhibit 
universal power-law scaling, with scaling laws suit ably 
modified by the presence of the trapping potential [l5|. 
With this in mind, we now turn our attention to prop- 
erties of the specific heat and compressibility in trapped 
Fermi gases in the vicinity of the superfiuid transition. 



III. SPECIFIC HEAT AND COMPRESSIBILITY 
IN TRAPPED GASES 

For a homogeneous three dimensional XY superfluid, 
thermodynamic quantities such as the specific heat Cy 
at constant volume and the compressibility kt at con- 
stant temperature exhibit simple scaling behavior near 



the critical point. For the specific heat, one has [11 1 



Cv = NkB [f\t\-'' + g + h,{n,T)] 



(8) 



Here, a ~ —0.01 [25|, |36|, and / and g are dimension- 
less coefficients, with different values above {fn,9n) and 
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below {fs,gs) the superfluid transition, and h{n,T) de- 
scribes the analytic background contributions to the spe- 
cific heat. In general, / and g depend on density but 
not temperature displays the asymptotic tempera- 
ture dependence in the limit T —^ T^). Similarly, for the 
compressibility kt = n^^{dn/d^)T, we have 



1 



ne-p{n) 



(9) 



where 7 ~ 1.3 [37] and k — k{n) is a dimensionless func- 
tion of density that assumes different values above and 
below the superfluid transition. hn{n, T) is a dimension- 
less function of density and temperature describing back- 
ground contributions to the compressibility. 

Equations ([8]) and ([9]) can be extended to inhomoge- 
neous systems by making use of the local density approx- 
imation (LDA): 



^i[n{r)] = /io - Vcxt(r). 



(10) 



Here, /i[n(r)] is the local chemical potential, /io is the 
equilibrium position-independent chemical potential, and 
V^xt is the external trapping potential. We can use pU]) 
to obtain an expression for the local superfluid transition 
temperature T'c(r) since T^. (x fic, the chemical potential 
at Tc. Thus, for a power-law potential Vcxtif) = Vpr^, 
Tc(r) = Tc(0)[l - xP], where x = v/Rtf and Rtf = 
{fiQ/Vp)p is the Thomas-Fermi radius of the gas. The 
local reduced temperature t{r) = [T ~ Tc{r)]/Tc{r) is 



tir) 



t{0) + xP 
1-xP ' 



(11) 



where t(0) = [r-rc(0)]/Tc(0) is the reduced temperature 
at the trap center. When t{Q) < 0, the system becomes 
superfluid, with the order parameter appearing at the 
trap center where t{r) is smallest. 

In the remainder of this Section, we apply the lo- 
cal density approximation (fTTj) to ^ and ([9]) to deter- 
mine the critical properties of the specific heat and com- 
pressibility in trapped gases. Such an LDA approach is 
standard in the cold atom literature (see, for instance, 
Refs. [l^lll, [3^), although we emphasize that it breaks 
down for very small t{0) when the coherence length be- 
comes too large. In Appendix |^ we give a brief discus- 
sion of the validity of LDA in the critical region. 



A. Specific heat in a trapped gas 

Using pT|) in ([8]) and recalling that the coefficients / 
and g assume different values above and below Tc, for 
small values of t{0), the specific heat per unit volume in 
a trapped gas is 



n{r)kF 



mr)\- 

fn\tir)\' 



'9s 

\-9n 



-hc{n,T) 
hh,{n,T) 



for 

for Rs 



r<R, 
<r<i?c. 



Here, Rs = i?TF|i(0)|p is the radius of the superfiuid 
region. When T < Tc(0), the system is superfluid for 
r < Rs and normal outside. Rc is the Ginzburg critical 
radius, defined by \t{Rc — Rs)\ = Gi. It describes the 
region is space where the local reduced temperature is 
inside the Ginzburg critical region and we expect the 
scaling behavior given by ([8]) to apply. In the normal 
phase, Rs — Q and the specific heat is given by the second 
line in p^ . 

Experiments 0, @ measure the bulk specific heat C = 
47r Jq^^" drr'^c{r) of a trapped gas. Using the above 
results and the Thomas-Fermi expression [l| n(r) = 
n(0)[l - for the density profile, for T < TdO) 

(i(0) < 0), we obtain the following expression for C: 

C = / dxx^f six, T)[l~xP]-^+''[m\-xP]-'' 

JO 

„2 J- rr^x n n]4+a 



+ / , dxx'Ux,T) [1 - .tP]'+" [xP - \t{ 
J\t(o)\p 

+H-{T). (13) 
For T > Tc(0) {t{Q) > 0), the bulk specific heat becomes 



C = 



K 



Ux, T) [1 - [xP + mr'^+H+{T). 

(14) 



(12) 



In these equations, R'^ = Rc/Rtf and / oc fR^p depends 
on n{x) and T (since -Rtf oc [/io(T)]p). Also, we have 
absorbed the contributions due to g as well as the back- 
ground contributions into the terms denoted by H^{T). 
These terms also include the noncritical background con- 
tributions which are not described by (fT2|) and come from 
the spatial region Rc < r < Rtf. 

Equations analogous to (fT3|) and (|14p have been used 
to determine the exponent a in high-precision measure- 
ments done on large sample sizes of superfluid ^He where 
gravity effects give rise to a position-dependent transi- 
tion temperature Tc(r) [I^IH, In these studies, the 
reduced temperature t(r) varied by only a small amount 
over the whole sample. This meant that one could be in 
a range of temperatures where the entire system was in 
the critical region and furthermore, in a single phase (i.e., 
superfluid or normal). Consequently, experiments could 
also extract information about the scaling amplitudes / 
and g from measurements of the bulk specific heat. In 
contrast, in trapped gases, t(r) varies significantly over 
the gas since the density decreases rapidly away from the 
center. Thus, the bulk specific heat of a trapped gas will 
always contain a contribution from the noncritical region 
at the edge of the gas. Furthermore, for t(0) < 0, one can 
never realize a region at finite temperatures where the en- 
tire gas is in the superfiuid phase. For these reasons, it 
is not possible to obtain detailed information about the 
scaling amplitudes / and g in trapped gases. 

Nonetheless, one can still use (IT^ and ([H]) to deduce 
the scaling behavior C oc (where a' is the exponent 
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for the trapped gas) of the specific heat in the critical 
region as long as the changes in H^{T) and Ho{T) are 
small over the range of absolute temperatures T probed. 
As with ^He, this requirement is easily realized since non- 
critical background contributions depend analytically on 
temperature. (The chemical potential of a unitary Fermi 
gas is also a smooth, slowly varying function of temper- 
ature [5]). What is different in inhomogenenous trapped 
gases is that the noncritical background terms H^{T) 
will almost certainly dominate over the critical contribu- 
tions to the bulk specific heat if the volume oc of the 
critical region is much smaller than the volume oc R^y of 
the gas. In this case, it would be difficult to experimen- 
tally resolve scaling behavior. A related issue arises in 
weakly interacting trapped Bose gases, where the contri- 
bution to the shift in (from its ideal gas value) due to 
noncritical background terms dominates over critical con- 
tributions as a result of the smallness of {Rc/RtyY 0- 

At T = 7c (0), the Ginzburg critical radius becomes 
Rc = i?TF [Gi I (Gi + 1)] ^ and the fraction of atoms in the 
critical region is thus ~ [Gi/{Gi + l)]p. For a Bose gas 
with scattering length a and thermal de Broglie wave- 
length At, above Tc, Gi - (a/Ar)^ < 1 [H. This means 
that, for atoms confined in a harmonic trap {p — 2), the 
fraction of Bose atoms that are in the critical region is 
of the order [40| {a/XrY ^ 1- contrast, for a Fermi 
gas at unitarity, the fact that Gi ~ C[l] means that 
Rc ^ OIRty] and a significant fraction of the gas is in 
the critical region when t(0) ^ 0. This suggests that scal- 
ing behavior of the bulk specific heat might be observable 
for strongly interacting Fermi gases. 

We start by considering the scaling behavior in the 
normal phase, T > Tc{0). Evaluating the integral in 
and expanding the result in powers of t(0), the leading 
terms are |4l| 



C^co + cii(0) + c2[t(0)]5-", 



(15) 



where cq, ci, and C2 are constants of the same order of 
magnitude. For the case of a harmonic trapping poten- 
tial {p = 2), the contribution to the specific heat com- 
ing from the third term in the above expression has also 
been obtained in Ref. [l^ (as one can see by applying 
the Josephson scaling relation [i^l a = 2 — 3i/ to the re- 
sults of Ref. [l^ , where v is the critical exponent for the 
coherence length). In addition, however, we find a term 
that depends linearly on the reduced temperature t(0). 
For harmonic traps, this term will dominate the specific 
heat of a trapped gas in the region about Tc(0). 

Below Tc(0), the situation is more complicated due to 
the inhomogeneity of the gas. In particular, note that 
the rate of change with respect to \t{0)\ of the integrals 
in the first and second lines of (|13p have opposite signs. 
Assuming that fs and /„ are both negative numbers (as 
is the case with '^He and superconductors), this means 
that while the slope of the specific heat still changes dis- 
continuously at T = T'c(O), the peak in the specific heat 
may occur below this cusp. Lowering the temperature 



below Tc(0), the specific heat may continue to increase 
until the first line in (|13p dominates over the second. The 
existence and size of the displacement of this peak from 
rc(0) depends on the values of the amplitudes /„, ^s, 
and gn ■ A similar situation arises in studies of inhomoge- 
neous ^He (see for instance. Fig. 2 in Ref. [1^). For these 
reasons, it is better to investigate the scaling behavior of 
the specific heat in the normal phase. 

Equation ([15]) shows how the scaling behavior of the 
specific heat of a uniform superfluid in the critical re- 
gion becomes strongly suppressed in a trapped gas due 
to the extra term 3 /p. (For positive critical exponents, 
the larger the exponent, the weaker the cusp.) It also 
shows how critical behavior is more pronounced in a 
"box-like" trapping potential: the higher the power of 
the potential, the smaller the critical exponent and the 
more the trapped gas behaves like a uniform superfiuid 
where C oc This suggests that one way to enhance 

critical effects is to trap the gas in as box-like a potential 
as possible. To this end, we also note that R'^ tends to 
unity as p oo, meaning that the fraction of atoms in 
the critical region is also larger in higher-order potentials. 



B. Compressibility in a trapped gas 

We now consider the behavior of the local compress- 
ibility KT{r) close to the critical point. The advantage of 
dealing with this quantity is that one can probe the local 
compressibility directly by looking at the density profile 
n{r) of the gas, instead of dealing with a trap averaged 
quantity such as the bulk specific heat. In principle, this 
allows for more direct access to the critical region. 

Using (fTO|) . one finds that the local compressibil- 
ity KT{r) = n~^{r){diJ,[n{r)]/dn)^^ of a harmoni- 
cally trapped (V^=2 = muj^/2) gas is related to the 
slope of the density profile by [sl, [3§] Krir) = 
— [mn?{r)uj'^r]~^dn(r)/dr. (The constancy of temper- 
ature is enforced by the fact that it is independent of 
position.) Combining this with ([9]) and (|lip . in the su- 
perfluid phase we find (restoring h) 



''ho 



dn{x) 



3 (a;) dx 



xn3 



m\ 



1 



(16) 



for the normal region r > R^ 
Here, 



{x > ^/\t{0)\) of the gas. 
afio = \/ h/muj is the harmonic oscillator length. 
A similar equation applies to the superfluid region r < 
Rs {x < ^|t(0)|). In (HI]), k oc {fio/huj)k and oc 
{fio/hu;)hi^ are dimensionless functions of n and T. 

Mean-field theories of trapped Fermi superfiuids (e.g., 
Ginzburg-Landau [1^ and BCS predict a disconti- 
nuity in the compressibility and hence, a "kink" in the 
density profile (where dn{x)/dx is discontinuous) at the 
local critical point x = yjt(0) using LDA [4j, |4l|. In 
contrast, the LDA scaling relation predicts that the 
slope of the density profile diverges at a; = yJUfi). Nei- 
ther of these features should be taken too seriously since 
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LDA breaks down in the immediate vicinity of the local 
superfluid transition. By looking at the region slightly 
away from the local critical point, however, detailed mea- 
surements of the density profile could in principle be used 
to determine the critical exponent 7. Wc recall that the 
fact that the Ginzburg critical region is large means there 
is a large spatial region where is expected to be valid. 

Without detailed knowledge of the background func- 
tion /ik(?t-, T) or the temperature dependence of the chem- 
ical potential fioiT), it is impossible to extract the crit- 
ical exponent 7 by measuring the density profile n{r) 
at a fixed temperature since Tp = Tp^n) is position- 
dependent and the value of the density n{r) can change 
rapidly over a small spatial region. However, by measur- 
ing the left-hand side of ([16]) for a fixed value of x close to 
the local critical point and varying t{0), for small changes 
in the absolute temperature T, it may be possible for ex- 
periments to resolve this critical exponent. (The weak 
temperature dependence of the chemical potential /io (T) 
close to Tc in the unitary region means that k is also 
weakly temperature dependent.) 



IV. IMPLICATIONS FOR EXPERIMENTS 

We turn now to the question posed in the Introduction 
and try and understand whether any signatures of critical 
behavior can be seen in the measured specific heat 0, 
and density profiles [i^ of a unitary Fermigas. 

The experiments reported in Refs. 0, 0| observe an 
anomaly in the specific heat at the putative transition 
temperature, but lack sufficient resolution for any strong 
conclusions to be drawn about the scaling behavior above 
Tc- It is thus important to estimate the range of tempera- 
tures above Tc where we expect the scaling law in (|15p to 
be valid and also the strength of the nonanalytic critical 
behavior (i.e., the size of ci and C2 in (jl5|) V 

The Ginzburg-Levanyuk number gives an estimate 
of the temperature range over which critical behavior 
manifests itself and we expect scaling behavior of the 
form given by ([5]) and ^ to be valid over much of 
this range. To check this point, we look to experi- 
mental results in other strongly interacting three dimen- 
sional XY superfluids, ''He 4g| and cuprate supercon- 
ductors, notably YBaCuO [47 1. Experiments on ''He have 
exquisite control over the temperature and are able to 
probe the critical region in a tiny region about Tc- The 
most extensive experiments have measured the specific 
heat over eight decades: 10~^° < \t\ < 10"^ These 
experiments confirm that scaling behavior of the form 
given by ([8]) is valid for at least \t\ < 10^^. In measure- 
ments of the specific heat of the cuprate superconductor 
YBaCuO, the same three dimensional XY scaling be- 
havior has been observed in the region \t\ < 10~^, larger 
than the Ginzburg-Levanyuk number Gi. An important 
reason why the specific heat in these strongly interacting 
three dimensional XY systems clearly exhibit the scal- 
ing behavior ([5]) over the entire Ginzburg critical region is 



that the coefficient / is large compared to g and he, and 
critical contributions to the specific heat dominate over 
noncritical background contributions (see, for instance, 
Ref. [13). 

Extrapolating these results to a strongly interacting 
Fermi gas close to unitary (which has a comparable crit- 
ical region), the scaling law in ([TS]) should be valid over 
a large temperature range, at least 10^ ^Tc ~ 0.02Tf. 
This is close to the resolution in current experiments 
which, for entropy S, we estimate to be around AT ~ 
{kBAE/SEF)TF O.OITf from the data in Ref. [i|. This 
suggests that these experiments are on the cusp of being 
able to resolve the linear scaling behavior predicted in 
([T5|) above Tc(0), assuming that background effects do 
not overwhelm contributions from critical behavior. 

The observability of the scaling predicted by (fT5|) is, of 
course, dependent on the size of the coefficients ci and 
C2. If these are small relative to the background contri- 
butions, then critical scaling may not be seen in exper- 
iments. However, as pointed out in Sec. IIIII the large 
size of the Ginzburg critical region means that a signif- 
icant fraction of the gas is in the critical region when 
T ^ Tc(0), in contrast to the situation in weakly interact- 
ing Bose gases. This means that as long as the amplitude 
/ in the scaling law ([5]) for the homogeneous gas is not 
too small compared to background terms, then the am- 
plitudes ci and C2 for the inhomogeneous trapped gas will 
not be too small either. We cannot say anything about 
the magnitude of / compared to hc{n,T) since the ratio 
of these quantities is nonuniversal and will be particular 
to the unitary Fermi gas. However, a small value of / 
(such that critical behavior is washed out by background 
effects) would be in sharp contrast to other strongly in- 
teracting three dimensional XY systems, where specific 
heat scaling is clearly evident. 

Finally, we recall that the peak in the specific heat 
may be situated below the transition temperature Tc(0). 
The scaling behavior predicted by ^5]) thus cannot be 
extended right down to the specific heat peak. Equa- 
tion (|15p is only valid for T greater than the superfluid 
transition temperature rc(0), where the specific heat ex- 
hibits a (likely) weak nonanalyticity. However, since 
the superfluid region grows quite fast in a trapped gas 

{Rs = -RTF|i(0)|i) and ~ /„ [H], from ([131), we ex- 
pect that the peak in the specific heat will be very close 
to T = Tc(0). Asymptotically, asp 00, the peak moves 
closer to T = Tc(0), giving rise to the expected nonana- 
lytic peak in the specific heat (lambda curve) character- 
istic of a uniform superfluid. 

Turning now to our prediction (|16p for the density pro- 
file in the vicinity of the critical point, we note that until 
very recently, no experiment on a trapped Fermi gas close 
to unitary has observed any feature in the density pro- 
file indicating the presence of a superfluid component. 
This suggests that the value of k is very small and the 
divergence in the compressibility is very weak. Note that 
this would be consistent with the situation in ^He, where 
the divergence of the isothermal compressibility is also 
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weak [H. 

Recent detailed studies of the density profile of a uni- 
tary Fermi gas carried out at MIT do reveal a feature 
in the density profile (see Fig. 48 in Rcf. ^45*|). How- 
ever, the resolution in these experiments is not yet good 
enough-and the effect is too small-to be able to measure 
the critical exponent 7. 

V. CONCLUSIONS 

In summary, we have analysed the width of the criti- 
cal region in a Fermi gas through the BCS-BEC crossover 
and discussed signatures of critical behavior in the spe- 
cific heat. We find that close to unitarity, the critical 
region is very wide-wider than clean cuprate supercon- 
ductors and comparable to superfluid ''He. Inside the 
Ginzburg critical region, any method for calculating ther- 
modynamic quantities that relies on an expansion in pow- 
ers of fluctuations of the order parameter is intrinsically 
unreliable and cannot be expected to produce quantita- 
tively accurate results. 

This is an important point to emphasize in light of re- 
cent comments (49j in the literature concerning unphys- 
ical features predicted by various Gaussian fluctuation 
theories, notably a double- valuedness of the order param- 
eter close to Tc d, m, [13] ■ This pathology is well-known 
in the Bose gas literature (see, for instance, Refs. f3TI.[5l|) 
and directly reflects the break down of perturbation the- 
ory in the critical region. Indeed, this double-valuedness 
occurs through the BCS-BEC crossover over a temper- 
ature range of the order Gi about (see Fig. 1(a) in 
Ref. [35J ) . One can thus interpret the numerical re- 
sults reported in Refs. [1, [35| as "empirical" evidence 
of a large critical region close to unitarity. Further ev- 
idence of a large critical region close to unitarity was 
found in Ref. [52| , where a diagrammatic analysis showed 
that diagrams that are subleading for T ^ Tc exceed 
leading ones in a wide region 0[Tc] above Tc- These re- 
sults emphasize that the only quantitatively reliable ana- 
lytic method in the critical region is the renormalization 
group p,[5i,[5i,[5|. 

Despite the fact that the width of the critical region 
of a unitary Fermi gas is of the order unity, signatures of 
critical behavior in trapped gases are suppressed by their 
inhomogeneity. For harmonic traps, we find that the bulk 
specific heat is linear in the reduced temperature close to 
Tc. Furthermore, the strength of the critical behavior is 
diminished in traps where the bulk specific heat includes 
contributions from the noncritical region at the edge of 
the gas. These factors likely explain why experiments 
have failed to observe a clear lambda-like curve in the 
specific heat, although a distinct feature is observed close 

toTe 0,i. 

We have also considered a scheme to measure critical 
behavior in the compressibility by looking at the den- 



sity profile. This scheme has the advantage that one can 
restrict the measurement to the critical region since it in- 
volves a local rather than bulk thermodynamic quantity. 
However, the nonanalytic contributions to the compress- 
ibility seem to be small in comparison with noncritical 
background contributions, meaning that experiments will 
need good resolution to be able to measure the critical ex- 
ponent for the compressibility. Nonetheless, a feature has 
already been observed in the density profile of a trap ped 
unitary Fermi gas below the superfluid transition [45| and 
future experiments may be able to explore this in greater 
detail. 
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APPENDIX A: VALIDITY OF LOCAL DENSITY 
APPROXIMATION 

Close to the superfluid transition, the local density ap- 
proximation (LDA) used to obtain the results in Sec. Hill 
should be valid as long as the length scale i?^ over which 
the local (Ginzburg-Landau) coherence length ^(r, T) 
varies significantly is much larger than the coherence 
length ^(r = 0, T) at the trap center: C(0,T)/i?j < 1. 
This is analogous to the usual finite-size criterion for a 
uniform superfluid confined in a box with sides of length 
L that requires £ /L <C 1 in order to be able to ignore 
finite-size effects (57| . For a gas confined in a harmonic 
potential V — the local coherence length is 

f(r,r) cx |r-Tc(r)|--'^ cx (1 - rVi?|)-'', where ~ 0.67, 

and = i?TFv|i(0)|- The Thomas-Fermi radius is 
(apart from an interaction renormalization factor of the 
order unity and restoring h) [l| Rtf = \/h/muj{24:Ny/^, 
where A'' is the atom number. Making use of the fact that 
the T = coherence length in a homogeneous unitary 
Fermi gas is [23l.[2^ ~ kp^, we take the coherence length 
at the trap center to be ^(r = G,T) - kp^{r = 0)|i(G)|"'', 
where kpir = 0) = y/2mEF/n'^ and Ep = hu}{my/^ is 
the bulk Fermi energy of a trapped gas. Thus, the condi- 
tion ^(0) < i?5 is satisfied when |t(0)| < (24iV)-2/(3+6<^), 
In typical experiments {N ~ 10^), this means that LDA 
is valid except in a small region \t{0)\ < 10-2 about i;. 
Note that this is equal to the estimate given in Ref. 
for a weakly interacting dilute Bose gas. For smaller val- 
ues of \t{0)\, when ^(0,r) ~ "finite-size" effects will 
become important. 
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